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Abstract 

The aim of the present paper is to introduce and to discuss inconsistencies 
errors that may arise when Eulerian and Lagrangian models are coupled 
for the simulations of turbulent poly-dispersed two-phase flows. In these 
hydrid models, two turbulence models are in fact implicitely used at the 
same time and it is essential to check that they are consistent, in spite of 
their apparent different formulations. This issue appears in particular in 
the case of very-small particles, or tracer-limit particles, and it is shown that 
coupling inconsistent turbulence models (Eulerian and Lagrangian) can result 
in non-physical results, notably for second-order fluid velocity moments. This 
problem is illustrated by some computations for fluid particles in a turbulent 
channel flow using several coupling strategies. 
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1. Introduction 

Polydispersed turbulent two-phase flows are found in numerous environ- 
mental and industrial processes, very often in conte xts that invo l ve ad ditional 



issues, for example chemical and combustion ones (jClift et al. 1 . 119781 ). From 
a practical point of view, the Navier-Stokes equations (for the fluid) and 
the particle equations (for the particulate phase) must be solved. Gener- 
ally speaking, two basic approaches have been used for the description of 
the particle phase: the Eulerian, or two-fluid approach, where the dispersed 
particle phase is treated as a fluid in much the same way as the carrier 
phase, namely by a set of continuum equations which represent the con- 
servation of statistical means, such as mass, momentum and energy, within 
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some elemental volume of the dispersed phase and the so-called Lagrangian, 
or particle tracking, approach where individual particles are tracked through 
the computed fluid field by solving the individual particle equation of mo- 
tion. In the latter case, the complete method for the two phases consti- 
tutes a hybrid Eulerian/Lagrangian approach. It is worth remembering 
that various strategies have been explored to couple different Lagrangian 
tracking approaches with Eulerian appro aches fo r the fluid phase, in par 



ticular, with RANS /Mom ents ap proach (IStock . 11994; iMuradoglu et al. 



200 it iMinier et al. 



LES 



Okong fc Bellan 
2006; Toschi fc Bodenshatz 



.. .2004) 

20041). and DNS 



2000; IJaberi et al. [ 2002; 



Soldati fc Marchioli I . 120091 ; iBec et al 



20091 ; lEaton 1 . 120091 ) . 



Let us now introduce the general nature of the Lagrangian approach. A 
no- model approach, in the spirit of DNS, is possible, but, in practice, the 
exact equations of motion are not treatable in realistic cases. Indeed, in the 
case of a large number of particles and of turbulent flows at high Reynolds 
numbers, the number of degrees of freedom turns out to be huge and one has 
to resort to a contracted probabilistic (modeled) description. In this case, 
particles are represented by an ensemble of Lagrangian stochastic particles 
whose properties are driven by either a model given i n the form of a set 



of stochastic differential equatio ns (continuous SDEs) (IMinier fc Peirano 



200ll ; IChibbaro fc Minier 1. 120081) or direct l y in terms of a nume r ical s cheme 



(random walk) (IStock I . I1994J : IPope 1 . 119871 : Maclnnes fc Bracco I . Il992h . The 



solution of the set of stochastic equations represents a Monte Carlo simula- 
tion of the underlying pdf. Therefore, this approach is equivalent to solving 
directly the corresponding equation for the pdf in the corresponding sample- 
space. In this work, we shall use the continuous approach. It is important 
to underline here that this choice is made for the sake of clarity and without 
loss of generality, since the Langevin approach is expressed in terms of con- 
tinuous variables and, thus, is physically more intuitive. The random walk 
models share the same properties but they are discrete and therefore they 
correspond to a numerical scheme for a given continuous stochastic model. 

In turbulent two-phase flows, the SDEs equations of the model contain 
several mean fields and have the general form 

dZi{t) = A(t, z, (/(Z))) dt + z > (/( z )» <Wj(t)> (l) 

3 

where the operator ( ) stands for the mathematical expectation. For exemple, 
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a typical Langevin model has the form (IMinier fc Peirano 1 . 120011 ) 



dx Pti = U v .idt (2) 
dU pA = -{U s>i - U Pti )dt (3) 

Tp 

dU ~ = -^Tr dt+i{u "^ {u " ))8J &f dt 



jL (U,,i ~ (Uf,i)) dt + ^ ( £ > (c„hklk + ^(k'k/k - 1) ) dW, (4) 



where quantities such as TJ i; k etc. are defined precisely elsewhere ( IMinier &: Peirano 



20011 ; IMinier et al. 1 . 120041 ). For the sake of present discussion, the important 
point is that this form reveals that different mean fields enter the model equa- 
tions. Fluid mean fields, typically (Ufj), are provided by the Eulerian solver 
while particle mean fields, such as (U p> i), are extracted directly from the par- 
ticles and are therefore provided by the Lagrangian solver. As such, in the 
case of particles with a non-negligible inertia, the problem is well-posed since 
the different mean fields come from different sources (Eulerian /Lagrangian 
solvers). In that situation, it may even be tempting to believe that improving 
the prediction of one mean field, for example the fluid mean velocity, results 
in improving the overall capacity of the complete model. 

However, the present hybrid method raises issues of consistency between 
the Eulerian and the Lagrangian solvers. This issue is particularly apprecia- 
ble in the tracer limit (fluid particle), for which the stochastic two-phase flow 
model reduces simply to a fluid model. This asymptotic limit-case represents 
the most relevant and fundamental situation where to test the effect of pos- 
sible errors due to inconsistency, for two main reasons: (1) in this limit, as 
will be shown below, possible inconsistency can be seen directly at the level 
of fluid moments; (2) one of the main application in bounded flows concerns 
particle deposition, where tracer particle are precisely the most important 
(aerosols) and, often, the least-well predicted by standard models. In the 
particle tracer limit of vanishing inertia, the model ([1]) takes the form: 

dx f>i = Ufjdt (5) 



P dxi Tj 



]r{Ui ~ (U t ) E )dt + y/C^jdWi , (6) 

L 



where Tl represents the Lagrangian time-scale of velocity correlations and 
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it is defined by T L = n /9 + 3M r7 , y 77\ ■ This model corresponds to the Simpli- 
fied Langevin Model (SLM) (IPope I . Il994bl ). In the above Langevin model, 
the fluid mean velocity field appearing in the rhs of the equation, (Ui) E , 
is provided by the Eulerian solver as indicated by the index E. However, 
since we are now dealing with fluid particles, the Lagrangian mean velocity 
field extracted directly from the particles, (Ui) L , represents the same physical 
property. We are therefore in presence of a duplicate field and the consistency 
issue requires that (Ui) L = (Ui) E . Yet, these two mean fields result in fact 
from two different sources: (Ui) E results from the turbulence model chosen 
in the Eulerian solver whereas (Ui) L results from the Langevin model. Thus, 
it is not obvious a priori to know what happens when the Eulerian physical 
model and the PDF one are not consistent. For instance, coupling Eulerian 
mean fields computed through DNS with a Lagrangian model wh i ch is consis- 



tent with a given RANS model (IPope I . Il994bl ; Muradodu et al. I . Il999l . 1200 lh 



may introduce a bias error. Unfortunately, this point has not yet received 
any attention and it has been too quickly believed that a "better" mean field 
(Ui) E fed into the simple Langevin model would automatically bring about 
a "better" model. This rou te has been already u s ed in many works, no t abb. 



a better model, inis rou te nas been already u s ed m many works, notably 
in particle deposi t ion cases flKroger fc Dressings I. I2000J : iMatida et al. I.l200(| 



Tian fc Ahmadi . 2007; Dehbi 



2008; 



Parker et al. 1 . 120081 : IZhang fc Chen 



20091 ). Even though in some cases and with some specific models this pro- 
cedure may turn out to be valid, it should be taken with care in general. 
This procedure is based upon the idea that the Eulerian and the Lagrangian 
parts of a hybrid method are completely independent, in absence of two-way 
coupling. It is our purpose to review critically this method. 

The question addressed in the present work is: what is the consequence 
on the particle simulations of using two inconsistent turbulence models (the 
turbulence model used in the Eulerian solver and the one corresponding to the 
Langevin model)? In order to develop a quantitative example, the purpose 
of this work is to investigate numerically this effect in a turbulent channel 
flow which is taken here as a relevant engineering case. 



2. Model consistency issue 

Following the discussion in the introduction, it is worth recalling that, 
in terms of Eulerian mean equations, the SLM model is equivalent to the 
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following turbulence model flPope I . Il994al fbl: iMinier fc Peirano LlioOlh : 

m) = 

ld(P) 

d(U t ) 

2 



dxi 

d(Ui) 
dt 



+ (Uj) 



d(Uj) t d(uiUj) 

dXj 



+ 



p dxi 



d(uiUj) d(uiUj) d{uiUjU k ) \d(Uj) 

' (U k ) a ^ + — t^z — = -(uiU k )-^- - {ujUk}- 



dt 



dxi 



dx k 



T, 



dx k 

■(uiUj) + C {e)Sij. 



dxi 



(7) 
(8) 

(9) 



Using the expression retained for Tl in Eq. dSJ), the transport equation for 
the second-order moments can be re-expressed as : 



d{uiUj) x d(u iUj ) 



dt 



dxi 



- -{UiUk}-^ 1 \UjUk) ^ 



dx k 

3 
-( 

2 



dxi 



v 1 + o^o) ( (UiUj) - -kSi, 



-Me). (10) 



This s hows that the SLM corresponds to a — e Rotta model (IPope 



1994bl ). We have retained this simple version, namely the SLM model, which 



is consistent with usual Reynolds-stress models as a kind of sound basis for 
the numerical investigations on particle deposition though it is clear that, 
at least for the prediction of fluid mean quantities, this leaves room for im- 
provement by using more complex Langevin ideas. 



3. Numerical Results 

A stationary turbulent channel-flow is solved first with the Eulerian method 
and, then, with the PDF one. The computations have been performed for the 
case of fully developed turbulent channel fl ow at Re = u T h/v = 395. The 



channel flow DNS results of Moser et al. (IMoser et al. I . Il999l ) have been 
taken for comparison as reference data. For the simulations, the mesh used 
in DNS (128 points along y) is also used here, which assures a good inter- 
polation of the mean fields for the Lagrangian stochastic particles. Quan- 
tities designed with the upper-script + are non-dimensionalised with wall 
parameters. Being a particle-mesh method, particles move in a mesh where, 
at every cell, the mean fields describing the fluid are known. Generally 
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speaking, the statistics extracted from the variables attached to the parti- 
cles, which are needed to compute the coefficients of system ([5])-(|I)J) are not 
calculated for each particle (this would cost too much CPU time) but are 
evaluated at each cell center following a given numerical scheme (averaging 
operator). The stochastic equations ar e solved through a consisten t 



first-order, unconditionally stable numerical scheme (jMinier et al. 1 . 120031 ). 
No-slip and impermeability conditions are satisfied by imposing boundary 
conditions on the stochastic particles. Symmetry conditions are imposed at 
the other boundary placed a the center of the channel y = h/2. Along the x 
direction, periodic conditions are imposed. The details of th e numerical ap 



proac h have been exhaustively described in a recent article (IPeirano et al. 



2006|). 



3.1. Consistency 

The first set of numerical experiments has been performed coupling the 
Eulerian Rotta Reynolds-stress model with the Langevin Simple model pre- 
viously described. Identical results are expected, since the equations for the 
moments of first and the second order are the same for both models and, 
thus, consistency is assured. In figures [Ik., the mean-velocity profiles ob- 
tained with the present hybrid configuration are shown. In figures [lb, the 
profiles of second-order moments, that is the Reynolds stress tensor, are also 
shown, for the present hybrid and for the Eulerian calculations. The Eulerian 
and Lagrangian profiles are in a quite good agreement. In conclusion, the 
mean fluid velocity moments derived from the present PDF method are in 
agreement with those computed in the Eulerian configuration and thus with 
the physical expected values. Moreover, the global hybrid method Eulerian 
(Rotta-model)/Lagrangian (SLM) is demonstrated to be consistent. It is 
worth emphasising that the ingredients that have been necessary to reach 
this objective are: 

(i) Consistent physical model. 

(ii) Consistent numerical scheme 

(iii) Accurate global numerical method , concerning also the exchange of 
information from Eulerian solver to Lagrangian one. 

3.2. Hybrid Consistency error 

We can now study the global error which is possibly introduced by us- 
ing a hybrid method not completely consistent. Generally speaking, hybrid 
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Eulerian/Lagrangian methods are affected by different kind of errors due 
to: (1) spatial discretisation; (2) time discretisation; (3) the use of a finite 
number of particle and per cell (statistical and bias errors). All these errors 
have been made negligible in the following simulations, in order to isolate 
the "hybrid consistency" error. In both methods, a time-step of 10~ 4 s and 
a spatial-step of 10 _4 5 have been used, which have been shown in numerical 
simulations to be sufficiently small for our purposes. 5 * 10 4 particles have 
been employed, which can be considered high enough in the present case. 
For the sake of clarity, the configuration discussed in the last section, where 
Eulerian mean variables consistent with Lagrangian ones were used, will be 
called the standard configuration in the following. 



First, let us consider the following configuration: an Eulerian DNS (IMoser et al. 



19991 ) is now coupled to the present Lagrangian method (this means that the 
mean fluid velocity (U) E provided to the SLM model is given by a DNS cal- 
culation). In figure [2ti, the mean velocity computed by the PDF method in 
this configuration, together with DNS original profiles as well as the mean 
velocity computed by PDF method in the consistent i2y-PDF configuration, 
are shown. The mean velocity obtained in the DNS/PDF configuration is 
in good agreement with the DNS one and, thus, is strongly different from 
the one obtained in the standard configuration. In fact, the exact profile is 
now recovered. In figure [2b, the Reynolds stress are shown for the same con- 
figurations. The profiles are now dramatically changed in comparison with 
the standard results. The (uv) and the (u 2 ) components are strongly over- 
predicted. On the contrary, the other diagonal components do not present 
appreciable changes. This behaviour can be explained. As diagonal (v 2 ) and 
(w 2 ) components are essentially dependent on Lagrangian time-scale and on 
diffusion coefficient, they are independent of the Eulerian model and there- 
fore not affected by the consistency error. On the contrary, the cross-shear 
stress (uv) depends not only on the turbulent kinetic energy and the La- 
grangian time-scale but also explicitely on the gradient of the mean velocity. 
It is therefore much more sensitive to the prediction of the Eulerian mean 
velocity. The present mean velocity profile attains a higher maximum than 
in the standard configuration and, thus, it is much steeper. The shape of (u 2 ) 
is in turn a direct consequence of this behavior. It is important to stress here 
that DNS/PDF Reynolds stress profiles are unphysical; notably, while the 
negative peak of the cross-shear stress should be of the order of 1, it turns 
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out to be (uv) ~ —3. 



Finally, we analyse the results obtained by coupling the present PDF 
model to a low-Reynolds number RANS model known for his good perfor- 



mance in boundary flows, the v2f model (IDurbin 1.119911). In the presen t 
calculations, we use a refined version of the model ([Laurence et al. 1 . 120051 ). 
In figure [3^, we show the mean velocity computed in this hybrid configura- 
tion together with the DNS and standard PDF results. As previously noted, 
the mean velocity given by the PDF method is approximately equal to the 
Eulerian mean one provided to the PDF solver, in this case computed by v2f 
solver. Furthermore, this result is also in good agreement with DNS result. In 
figure [3b, the Reynolds stress profiles for the same configurations are shown. 
The results are similar to those obtained in the hybrid DNS /PDF configura- 
tion. Even though the (uv) and the (u 2 ) components are less overpredicted, 
they still show qualitatively the same behaviour and remain unphysical. The 
consistency error is slightly less important but still large. 
Some observations are in order: 

(i) In all configurations, the mean velocity computed from stochastic par- 
ticles basically collapses on the value given by the Eulerian mean ve- 
locity used in the Lagrangian model. This is in line with the physics 
of the Langevin model, wh ich is based on a return-to-equilibrium idea 



(IMinier fc Pozorski 1 . 119971 ). 



(ii) In the PDF-DNS and PDF-V2F configurations, a large difference be- 
tween the Eulerian and the Lagrangian results at the level of the second- 
order statistical moments (Reynolds stress) is found. This is a direct 
consequence of coupling, in a hybrid approach, two methods which are 
not consistent. Thus, the hybrid consistency errors are identified by 
comparing the actual Reynolds stress profiles with those computed in 
the standard configuration, see figs. [2bj3b. 

(iii) Results obtained in the hybrid DNS/PDF and v2f/PDF approaches 
are qualitatively similar, even though the DNS and the v2f approaches 
are quite different from a theoretical point of view. This should be ex- 
pected. The profiles provided by the Eulerian solver to the Lagrangian 
one are the mean velocity, the mean pressure, the turbulent energy and 
the turbulent dissipation. For these variables, the v2f approach gives 
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results which are in good agreement with the DNS ones. Therefore, 
from the point of view of the hybrid method DNS and v2f approaches 
are qualitatively similar. Moreover, from a quantitative point of view, 
the results obtained through the hybrid v2f/PDF approach are in bet- 
ter agreement with those calculated in the standard Rij/PDF case than 
those obtained using the DNS/PDF approach, fig. [3b. Therefore, v2f 
model is found to be more consistent with present Lagrangian model 
than DNS and the hybrid consistency error is smaller. 



4. Conclusions 

A study of a turbulent channel flow has been carried out in the framework 
of a Hybrid Eulerian/Lagrangian approach with the main attention pointed 
to the issue of the consistency between the two approaches. 

We have used different kind of Eulerian methods (RANS Rij and v2f, 
DNS) coupled with the simplified Langevin model for the Lagrangian tracking 
of fluid particles. It is worth underlying that this kind of couplings are nor- 
mally used in mul t iphas e simulat i ons, f or inst ance DNS (Matida et al. . 2000l : 



Tian fc Ahmadi I . 120071 : iDehbi I . 120081 ). v2f flZhang fc Chen I . 120091 ) among 



others. The Lagrangian model proposed is in the form of a set of stochas- 
tic differential equations ©-([6]) and, in practice, the PDF underlying this 
stochastic process is obtained via a Monte Carlo method, that is through 
the simulation of a number of stochastic particles. This feature is very use- 
ful to recognize immediately the physical contain of the model, at variance 
with other more heuristic but also very popular approaches like the discrete 
random walk models. 

We have analysed the changes produced in the Lagrangian results by 
choosing different Eulerian methods. In particular, we have studied the du- 
plicated variables which are the two first moments and we have shown that 
the global method is found to be consistent when Rotta model is coupled 
with the SL model, as it is theoretically expected. On the contrary, results 
obtained in other configurations give evidence of a bias error precisely due to 
the inconsistency between the Eulerian and Lagrangian physical models. It is 
worth emphasizing here that it is possible to assess the global consistency of 
the method, because we use a completely consistent and an accurate numer- 
ical scheme, which allows us to consider the numerical errors as negligible. 

These results underline that, in hybrid Eulerian/Lagrangian approaches 
to turbulent flows, delicate questions of consistency between the two parts 
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of the method ineluctably arise and deserve a careful treatment. In general, 
the simple statement that a "better fluid profile" help to improve Lagrangian 
results (independently of the Lagrangian model) is not true. At variance 
with this belief, using an Eulerian model which is not consistent with the 
Lagrangian one introduces a consistency error which can introduce a flaw in 
the global method. 

Finally, in authors' opinion, the Lagrangian tracking method is too often 
considered as a black box tool scarcely important from a physical point of 
view and thus hierarchically subjected to the Eulerian method. In fact, 
we think that in hybrid Eulerian/Lagrangian approaches to two-phase flows 
the Lagrangian part is more fundamental and has the guiding role. In this 
sense, our results show that it could be wise to start from the choice of 
the Lagrangian model and afterwards to choose an Eulerian model which is 
consistent with it. 
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